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I. INTRODUCTION 



The classical linear stability theory of viscous incompressible flows is concerned with 
the development in space and time of infinitesimal perturbations around a given basic flow. 
Then small disturbances are resolved into normal modes which, for a steady-state basic flow, 
depend on time exponentially with a complex exponent A. For parallel shear basic flows, 
further separation of variables in the governing stability equations leads to a set of ordinary 
differential equations which, with taking recourse to Squire's theorem and considering only 
2-D disturbances, reduces to the Orr-Sommerfeld equation. When this equation is solved 
with proper boundary conditions, the problem of linear stability of parallel flows is reduced 
to a 2-point boundary (eigen) value problem. 

For nonparallel basic flows, when the equations for disturbance flow are dependent not 
only on the normal coordinate, the corresponding operator does not separate unless certain 
terms are ignored. The approximation, that neglects the nonparallel terms and relates the 
stability characteristics to those of the equivalent parallel flow, has been extensively used 
for the boundary layer type flows to retain the great advantage of reducing the disturbance 
equations to ordinary differential equations (see, e.g., Reed and Saric, 1996). A number 
of weakly nonparallel theories, which seek to account for the affects of the flow divergence 
through equations that include higher-order terms than those in the Orr-Sommerfeld equa- 
tion, have been developed (see reviews in Reed and Saric, 1996; Herbert 1997; Saric et al. 
2003). 

It is worth mentioning, in this context, the works on stability of conical flows by Shtern 
and Hussain (1998), Shtern and Drazin (2000) where an exact transformation, reducing the 
stability problem to an ordinary differential equation despite the non-parallelism of the basic 
flow, is found. However, the advantage of such a transformation in the analysis is limited 
to the particular class of steady perturbation modes. In Shtern and Hussain (2003), the 
approach of Shtern and Hussain (1998) has been applied to the time-oscillatory disturbances 
by using a far-field approximation in the equations for the disturbances. 

Several methods have been designed for numerical solution of the nonparallel flow stability 
problems. In the context of boundary-layer type flows, the most successful effort to-date 
is the parabolic-stability-equation (PSE) approach, introduced and recently reviewed by 
Herbert (1997). However, the PSE based numerical studies are not able to accommodate 
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the upstream propagation of disturbances. The 'global' linear stability analysis (see review 
by Theofilis 2003) was developed for analyzing stability of the two-dimensional flows. The 
nonparallel stability effects have been also investigated on the basis of direct numerical 
simulations of the linearized Navier-Stokes equations (e.g., Davies and Carpentier 2003) or 
the complete Navier-Stokes equations (e.g., Fazel and Konzelmann 1990). Such numerical 
simulations are not equivalent to a stability analysis, and, in fact, have more in common with 
physical experiments than stability theory. The numerical methods suffer from the problem 
of boundary conditions on the 'open', inflow and outflow, boundaries which (especially 
outflow boundary conditions) can lead to spiritus effects, even when carefully implemented. 

If the basic flow is non-steady, this brings about great difficulties in theoretical studies 
of the instability since the method of normal modes in its traditional form, with the modes 
depending on time as exp(Xt), is not applicable at all. (Some success has been achieved in 
studying stability of the time-periodic basic states when Floquet theory can be applied - 
see Drazin and Reid 1981.) If an unsteady flow is non-parallel, it should further complicate 
matters. As a matter of fact, there are no examples of the linear stability problem for 
viscous incompressible flows developing both in space and time which is exactly solvable via 
separation of variables. 

Recently, the so-called direct approach to separation of variables in linear PDEs has been 
developed by a proper formalizing the features of the notion of separation of variables (see, 
e.g., Zhdanov and Zhalij 1999a, 1999b). In this approach, a form of the 'ansatz' for a solution 
with separated variables in a new coordinate system as well as a form of reduced ODEs, that 
should be obtained as a result of the variable separation, are postulated from the beginning. 
The method has been successfully applied to several equations of mathematical physics (see, 
e.g., Zhalij, 1999; Zhdanov and Zhalij, 1999a, 1999b; Zhalij, 2002). 

In the present paper, we apply this approach to the linear stability equations which govern 
the disturbance behavior in viscous incompressible fluid flows. The calculations are made 
for the linear stability equations written in cylindrical coordinates. The new coordinate 
systems and the forms of basic flows, which permit the postulated form of separation of 
variables in the equations for disturbances, are determined as the result of application of 
the method. Then the basic flows are specified by the requirement that they exactly satisfy 
the Navier-Stokes equations. 

The paper is organized as follows. In Section II, we give a description of the method, 
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as applied to the problem of linear stability of a three-dimensional unsteady basic flow 
with respect to the three-dimensional unsteady perturbations, and present the results. An 
example of application of the method to the linear stability equations with a restriction 
to the two-dimensional perturbations is also presented. In Section III, we discuss the fluid 
dynamics interpretation of some basic flows, defined in Section II as the exact solutions of 
the Navier-Stokes equations possessing exactly solvable stability problems, and formulate 
the corresponding two-point boundary value problems of ordinary differential equations. 
Concluding remarks on the results obtained are furnished in Section IV. 



II. APPLICATION OF THE DIRECT METHOD TO SEPARATION OF VARI- 
ABLES IN THE STABILITY EQUATIONS 

A. Procedure 

We formulate the linear stability problem based on the Navier-Stokes equations written 
in cylindrical coordinates (r, tp, z). As usual in stability analysis, we split the velocity and 
pressure fields (v r , v v , v z ,p) into two problems: the basic flow problem (V r ,V v ,V z , P) and a 
perturbation one (v r ,v^,v z ,p), 

v r = V r + v r , v ip = V ip + v v , v z = V z + v z , p = P + p (1) 

Introducing into the Navier-Stokes equation written in terms of the variables (v r , v^, v z ,p) 
and neglecting all terms that involve the square of the perturbation amplitude while imposing 
the requirement that the basic flow variables (V r ,V<p,V z , P) themselves satisfy the Navier- 
Stokes equations, one arrives at the following set of linear stability equations in cylindrical 
coordinates: 
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dv z dv z dV z V v dv z v v dV z dv z dV x 
+ V r — — h v r — h — -t; 1- — 1- V z — — I- v z - 
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Let us introduce a new coordinate system £, £ = r), 7 = 7(t, <f), rj — rj(t, z). 
We will say that the system (J2J) is separable in the non- stationary cylindrical coordinate 
system £, 7, 77 if the separation ansatz, 

f r = T(t) exp(ar/ + m-7 + sS(t))f(£), 

— T(t) exp(a?7 + mj + sS(t))g(£), (3) 
t> 2 = T(t) exp(a?7 + 7717 + sS(t))h(£), 
p = Ti(t) exp(a?7 + 7717 + sS'(t))7r(^) 

reduces the system of PDEs (J2J) to a system of three second-order and one first order ordinary 
differential equations for four functions f(£),g(£), h(£), 7r(£) of the following form 

h"(0 = U ll9 '(0 + U l2 h'(0 + U 13 tt'(0 + U 1A f(0 + U 1B g(£) + U 16 h(0 + U 17 ir(£), 
/"(£) = U 21 g'(0 + U 22 h\0 + U 23 tt'(0 + Uuf(0 + u 25 g{0 + U 26 h(0 + U 27 tt(0, (4) 
<?"(£) = U 31 g'(0 + U 32 ti(£) + U 33 ir'(£) + EWtf) + U x g(0 + U 36 h(£) + U 37 n(0, 
/'(O = ^4i/(0 + UvgiO + U i3 h(® + U u ir{S). 

Here C/y are second order polynomials with respect to spectral parameters a, s, m with 
coefficients, which are some smooth functions on £. 

Note, that equations ©-© form the input data of the method. We can change these 
conditions and thereby modify the definition of separation of variables. For instance, we 
can change the order of the reduced equations (jlj) or the number of essential parameters 
a, s, m. So, our claim of obtaining the complete description of basic flows and non-stationary 
coordinate systems providing separation of variables in (J2J) makes sense only within the 
framework of the given definition. If one uses a more general definition, it might be possible 
to construct new coordinate systems and basic flows providing separability of the system 

The principal steps of the procedure of variable separation in the system (J2J) are as follows. 
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1. We insert the ansatz (J3J) into equation (J2J) and express the derivatives /"(£)> 9"(£)> 
h"(£), f'(£), in terms of functions #'(£), tt'(0, /(0> MO; ^(0 usin S Ra- 
tions (gj). 

2. We regard g'(£), /i'(0> ^'(O? /(0> #(0; M0> n (0 as ^ ne new independent variables. 

As the functions £(t, r), j(t, <p), r)(t, z),T(t),Ti(t), S(t), basic flows V r , V 9 , V z and coeffi- 
cients of the polynomials Uij (which are some smooth functions on £ ) are independent 
on these variables, we can demand that the obtained equality is transformed into 
identity under arbitrary g'(£), h'(£), /(C)' MO? MO? MO- I 11 other words, we 

should split the equality with respect to these variables. After splitting we get an 
overdetermined system of nonlinear partial differential equations for unknown func- 
tions r),j(t, ip),r](t, z),T(t),Ti(t), S(t), basic flows V r , V v , V z and coefficients of the 
polynomials U^. 

3. After solving the above system we get an exhaustive description of coordinate systems 
providing separability of equations (J2J) according to our definition. 

Thus, the problem of variable separation in equation (J2J) reduces to integrating the overde- 
termined system of PDEs for unknown functions £(£, r),"f(t, <p), r)(t, z),T(t),Ti(t), S(t), basic 
flows V r ,V v , V z and coefficients of the polynomials U^. This have been done with the aid of 
Mathematica package. 

B. Results 

We will consider the stability problems with separated variables for the basic flows spec- 
ified by the requirement that they exactly satisfy the Navier-Stokes equations. 

1. Three-dimensional perturbations 

The forms of the perturbations v r , v v , v z and p are: 



v r = T(t) exp (^ar] + rrup + s J T(t) 2 dtj /(£), 
v v = T{t) exp ^ar] + rrup + s J T(t) 2 dtj g(£), 
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v z = T(t) exp ^ar] + mip + s J T{t) 2 dt^j h(£), 

where 

Z = T(t)r, V = T(t)z + c(t). (6) 

Two classes of basic flows satisfying the Navier-Stokes equations are found as the result 
of the analysis. Velocity fields for both classes are defined by 

T'(t) 

Vr = B^)T(f)-r^, 

V v = C(OT(t), (7) 
where the functions T(t) and B(£) are specified in different ways for the two classes. 

Class I: T (t) = ^, S(0 = -^ + |, (8) 
where the functions A(£) and C(£) satisfy the equations 

(4k + 3£ 2 - 4v)A(£) + f (-4A; + 3£ 2 + 4u)A"(£) + 4£W(0 = 0, 
-4z/fc £ + (-4A; + 3£ 2 - 4i/)C(£) + f (-4A: + 3£ 2 + 4i/)C"(0 + 4i/£ 2 C"(0 = 0. (9) 

and the pressure distribution is given by 



P vk ip x 



2 



p t 8t 2 



+ 1J W 3 dt + Po® ( 10 ) 



The ODEs Q can be explicitly solved in terms of the incomplete gamma functions. 

Class II : T(t) = 1, = - 

where A(£) and C(£) satisfy the equations 

(A; - ^)A'(0 + £(v - k)A"(0 + £W(£) = 0, 
vk £ + ( fc + ^(7(0 + £(* - i/)C"(0 - £VC"(£) = (11) 



and the corresponding pressure distribution is 



P 



P 



vko<p + x[f3\t) + vA!'{£) + 



v — k 




The ODEs can be explicitly solved in elementary functions. 

The equations with separated variables can be written for both classes in the forms 



2(mv - £C(0)g(0 + £((-!/ + £B(0)/'(0 + t{At) ~ "AO)) = 0, 

f(0(-2mu + e<7(o + mo) + e(m7r(e) + (-^ + womo - £*v'(£)) = o, 

/(0 + mp(0 + eWO + /'(0) = o- 



£ Two-dimensional perturbations 

The stability properties of a given flow may be tested by considering perturbations of 
specific structures. For example, the problem may be restricted to the two-dimensional 
perturbations even though an analog of the Squire theorem cannot be proved (see, e.g., 
Griffond and Casalis 2001, Joslin 1996), or the perturbation flow field may be taken to 
have the same general form as the basic state (Duck and Dry 2001). Although the stability 
analysis restricted to perturbations of specific forms is not complete, it enables one to show 
that the flow is susceptible to a special kind of instability. To demonstrate that a specification 
of the disturbance field may lead to new possibilities we consider the results of application of 
the direct method to the linear stability equations with a restriction to the two-dimensional 
perturbations with v z — and v r and v$ not dependent on z. 

The separability analysis leads to the perturbations of the form 




v z = 
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p = p T(t) 2 exp (^rmp + s j T(t) 2 d?j vr(£), 

Z = T(t)r (12) 

which is a particular case of (J3J) for a = 0. However, for the perturbations of the form ()12|). 
the corresponding basic flows are not restricted to those listed in Section IIB1. In addition, 
the following basic flows are permitted 

V z = -kz + /3(t), 
V r = kr/2 + q/r, 
V 9 = vB(£)T(t). 



- p = -^v + * (w - m)) - + nt) (V / ^±dt - \ 



where the functions T(t) and B(£) satisfy the equations 

T'(t) - i (gr 3 (t) - kT(t)) = o (13) 

k £ - (2g + 2z/ + ge 2 )^(0 - £(2<? - 2z/ + Qf 2 )B'(0 + 2^ 2 J B"(£) = 
which leads to the following cases 



Ve** + 1 / Ve fct - 1 V fc 



T(t) = l (| = l), T{t) = e- kt / 2 (Q = 0). 



III. SPECIFIC FLOW STABILITY PROBLEMS 



In this section, we will discuss the fluid dynamics interpretation of some basic flows, de- 
fined above as the exact solutions of the Navier-Stokes equations for which the correspond- 
ing stability problems are exactly separable, and will formulate the corresponding two-point 
boundary value problems. 

We will consider particular cases of the class of the exact solutions of the Navier-Stokes 
equations in cylindrical coordinates identified in Section IIB1 as Class I. It is possible to 
enrich the solution defined by the formulas © - (fTUj) using the invariance of the solution 
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with respect to a shift of time variable. Making change of variables t = t' — 1/6, where b is a 
constant, and omitting primes in what follows, we will have the solution of the Navier-Stokes 
equations in the form 

bz F(C) _ 1 
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+ -SQr=-/9(f), V r 
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where b can be both positive and negative. 
The corresponding pressure distribution is 
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The functions F(() and M(() satisfy the equations 



(15) 



(4* _ 3b( 2 - 4v)F'(() + C(-4* - 36C 2 + 4i/)F"(C) + 4CVF'"(C) = 0, (16) 
4Jfe„i>C + (-4Jfe - 3&C 2 - 4z/)M(C) + CM* - 3&C 2 + 4z/)M'(C) + 4CVM"(C) = 0. (17) 

Equations (fTHjl and (fTTjl can be solved in quadratures 

F(C) = ci + c 2 r (J-, -Z(C)) + c 3 y e^C*"^ (l - 1;, Z(C)) d£ 



Z(C) = 
M(C) 



3&C 2 



57/ 



(18) 



c 4 + c 5 r ( i + — , -z(c) 



fco f 3& 
2 V8f 



2 i.' 



where r(A, Z) is the incomplete Gamma function and c\, 
Note also the expression for F'(Q 



(19) 

C5 are arbitrary constants. 



F'(C) = 



-2c 2 (-!)%c 3 r v i--,z(c) 



(20) 



The correspondingly specified perturbations © take the forms 

v r = (1 — bt) s exp (arj + rrup) /(C), 
u v = (1 — fr£) s exp (ar/ + m<£>) 
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p 

where 



(1 — bt) s exp (arj + rrup) h((), 

(l-bt) s - 1/2 exp(ari + m V )7r((). (21) 



" = 7ra +c(() (22) 



The equations for the perturbation amplitudes are 

(-4k + b( 2 - 4b( 2 s +4v- Am 2 u - 4a 2 ( 2 v + 4a( 2 F{() + 4m(M{()) /(C) 

+C (4* - 4i/ + 3b( 2 ) f'(C) - 4K 2 /"(C) 

+8 (rm/-CM(C))y(C) + 4CV(C) = 0, 
(4fc + 6C 2 - 46C 2 s + 4i/ - 4m 2 z/ - 4a 2 C, 2 v + 4a( 2 F(C) + 4mCM(C)) ^(C) 

+C (4A: - 4z/ + 36C 2 ) ^(C) - 4i<Y(C) 

+ (_ 8mzy + 4C M(C) + 4C 2 M'(C)) /(C) + 4mC7r(C) = 0, 
-2 (6C 2 + 2&C 2 s + 2m 2 z/ + 2a 2 C 2 ^ - 2a( 2 F{() - 2mCM(C)) h(() 

+C ( 4fc _ 4z/ + 36C 2) h / (c) _ 4z/C 2 ^(C) + 4C 2 F'(C)/(C) + 4aC 2 7r(C) = 0, 
/(C) + mg(C) + C (ah(C) + /'(C)) = 0. (23) 

This system can be reduced to a system of two third-order equations for two functions (/(C) 
and g{(), for example). 

The above formulas (|14 p -([23| ) remain valid if we introduce the nondimensional variables, 
with the time scale l/\b\ and the correspondingly defined velocity scale. In the dimensionless 
equations (we will retain the same notation for the nondimensional variables), the parameter 
b takes one of the two values: b — 1 or b — — 1, and v is replaced by 1/Re where Re is the 
Reynolds number. 

Note that the solution (|T4^) - (JT5j) for b = 1 undergos the finite-time 'breakdown' (see, e.g., 
Duck and Dry 2001, Hall et al. 1992) and the 'normal mode' forms (j21j) -(|2*2* j) are naturally 
adjusted to the description of the disturbed flow as the breakdown time t — 1 is approached. 

Next we will consider several two-point boundary value problems corresponding different 
specifications of the basic flow (JHJ)-([Tnj). 

(i) Axially symmetrical stagnation-point type flows. These are the simplest basic 
flows, that are obtained from ()14|1 - ()15j) by setting 

F(C) = 0, M(C) = 0, c(t) = 0, k = 0. (24) 
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FIG. 1: Unsteady axially symmetric stagnation point flows: 6 = 1 (left), b = —1 (right) 

In the case of b = 1, the solution describes impingement of two axially opposite stagnation 
point flows with velocities growing with time, and, in the case of b = —1, the solution 
describes the flow where fluid flowing radially from infinity approaches the axis and spreads 
along it, with the flow velocity decreasing with time as (1 + t) _1 (Fig. Q). In both cases, the 
boundary conditions to equations (J2*3*)) are set at the axis ( = and at ( = oo, as follows 

/(0) = 0, g(0) = 0, h'(0) = 0; /(oo) = 0, g(oo) = 0, h(oo) = (25) 

(ii) Flow outside an expanding cylinder. This case corresponds to b = — 1. The radius 
of the cylinder changes with time as R = a/1 + 1. (In the dimensional variables, marked 
with stars, it is R* = Rq\/1 — b*t* and the value Rq is used as a length scale for the 
nondimensional variables while the time scale is 1/|6*|). The surface of the cylinder stretches 
in the longitudinal direction according to the law U = Kz where K = ~(1 + 1)^ 1 and U is 
an axial velocity at r = R. 

If the cylinder surface is impermeable, then using the boundary condition V r — V at 
r = R, where V = dR/dt = |(1 + 1) -1 / 2 is the radial velocity of the cylinder surface, yields 
k = 3/4. In the case of porous cylinder, we have k = 3/4 — Vo, where Vq is a constant defining 
a magnitude of the suction (Vo > 0) or injection (Vo < 0) velocity as V s = V (l + 1) -1 / 2 

Since the domain is infinite in radial direction, it should be set C3 = in the expression 
(I2(J|) for -F'(C) riot to have an unbounded behavior for F(Q at infinity. (The incomplete 
Gamma function T(A, Z) ~ Z A ~ 1 e~ z as Z — » 00 which results in F'(Q ~ as £ — > 00 and 
produces the logarithmic behavior for F(() at infinity - this can be confirmed by considering 
a behavior of F(() itself for specific values of k = 2nu, with n being a positive integer, when 
closed- form solutions of equation (fTBj) for F(Q can be found.) Then the flow at infinity 
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represents a combination of a stagnation point flow and a uniform stream, and the two 
constants in the expression for F(() are determined from the boundary condition V z — U at 
r = R {F(l) = 0) and the condition for the uniform part of the flow velocity at infinity be 



of the form U oc /\^l + t where is a constant. 

For not swirling flows (a swirl can be also added with the swirl velocity defined by (j!4)) and 
(HHJ)) we have the following to be introduced into the equations for perturbations (|2*T|) - (j2^j) : 



-n / fcRe 3Re( 2 
1 2 ' 8 



-p | fcRe 3Re 
1 2 ' 8 



M(C)=0, c(t)=0, k = ^-V Q . (26) 



where Re = _Rq 2 |6*|/z/* is the Reynolds number (the corresponding flow structure is illus- 
trated by Fig. EJ). Note that the solution (|2*Bjl is expressed in elementary functions for 
Re = 2n/k, where n is an integer, with the use of the specific values of the incomplete 
Gamma function (Abramowitz and Stegun 1965): 

r M = i-( 1 + z + - + ... + — y* 

In the equations for perturbations (}2~T]) - (f2lH) . it should be also set b = — 1 and v = 1/Re. 
The boundary conditions to equations ()23|) are set at £ = 1 and at ( = oo, as follows 



/(l) = 0, g{\) = 0, h{\) = ; /(oo) = 0, g(oo) = 0, h(oo) = (27) 
( iii) Flow inside an expanding porous cylinder. In this case, like as in the previous 



one, b = — 1, the radius of the cylinder changes with time as R = a/1 + 1 and the surface 
stretches according to the law U = Kz where K = 1(1 + 1) -1 . The difference is in that the 
fluid is now inside the cylinder and there is an injection of fluid through the porous pipe 
surface, which may be either normal to the surface or oblique, with the blowing velocity 
varying with time as = Vq(1 + 1)" 1 ^ 2 , where Vq is a constant. 

The conditions V r = at r = requires k = and using the condition dV z /dr = at 
r = in flUD yields 

3C 2 Re 

where Re is the Reynolds number. Then two arbitrary constants in the expression for F(Q 
are determined from the condition at the cylinder surface F(l) = (3/4)tanh#, where the 



13 




14 



angle 9 defines the direction of blowing (with respect to the inward radial direction), and 
from the condition at the axis F(0) = Uq, where Uq is a constant defining the axial flow 
velocity. Restricting ourselves to not swirling flows and normal blowing (8 = 0), we have 
the following to be introduced into the equations for perturbations (|21)) - (J23j) : 

Ei -Ei(-fs) _i nC 2 
HO = U \ + r( o ; 3| g)+ln( 3| g) > M(C) = 0, c(t) = 0, k = 0. (28) 

where Ei(Z) is the exponential integral function and 7 is Euler's constant (the corresponding 
flow is shown in Fig. |3J). Note that, despite the presence of the logarithmic term in the 
nominator, the expression (J28|) for F(() is finite at ( = since the expansion of Ei(Z) for 
small Z includes the term In Z. 

It should be also set b = — 1 and v = 1/Re in equations (|2*Tj) - (|2l?|) . The boundary- 
conditions for the perturbation amplitudes are set at the axis ( = and at the cylinder 
surface ( = 1, as follows 

f(0) = 0, g(0) = 0, h'(0) = ; /(l) = 0, g(l) = 0, h(l) = (29) 
(iv) Flow inside a contracting porous cylinder. In this case 6 = 1 and the radius of 



the cylindrical tube changes with time as R = a/1 — t (t < 1). The surface of the tube shrinks 
according to the law U = —Kz where K = |(1 — t)^ 1 and there is a suction of fluid through 
the permeable cylinder surface, which may be either normal to the surface (F(l) = 0) or 
oblique (F(l) = (3/4)tanh# with the angle 9 defining the direction of suction with respect 
to the outward radial direction). The suction velocity varies with time as VJ, = Vq(1 — t)^ 1 ^ 2 
where V$ = [9/16 + F^l)] 1 ^ 2 . For not swirling flows and normal suction (9 = 0), we have 
the following to be introduced into the equations for perturbations (PT|) - (f2^|) : 



Ei (K^.Ei^-lnC 2 



nO = U ^ _ Ei (3|.) + In (f^) ' M ^ = °' C ^ = °' k = Q - (3°) 

where Uq is a constant defining the axial flow velocity (the corresponding basic flow is shown 
in Fig. EJ). 

It should be set b = 1 and v = 1/Re in equations (|2*T]) - (f2l?j) . and the boundary conditions 
are set at the axis £ = and at the cylinder surface ( = 1, as follows 



f(0) = 0, g(0) = 0, h'(0) = ; f(l) = 0, g(l) = 0, h(l) = (31) 
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c) 




FIG. 3: Flow inside an expanding porous cylinder for Re = 100 and U = 5 at different time 
moments: a) t = 0; b) t = 1; c) t = 5. 

(v) Flow in the gap between concentric cylinders. Here, like as in the previous 
cases, different boundary conditions can be considered. 



IV. CONCLUDING REMARKS 

Several classes of the exact solutions of the Navier-Stokes equations describing spatially 
developing and unsteady flows, for which the linear stability problems can be rigorously 
reduced to eigenvalue problems of ordinary differential equations, have been defined. Those 
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FIG. 4: Flow inside a contracting permeable cylinder for Re = 100 and U = 6 at different time 
moments: a) t = 0; b) t = 0.5; c) t = 0.75. 

exactly solvable nonparallel and unsteady flow stability problems can provide a necessary 
foundation for a number approximate approaches used in the stability analysis so far. The 
results can be also used for testing the methods based on direct numerical simulations 
of the (linearized) Navier-Stokes equations. Note that the basic flows considered in the 
paper belong to the category of the so-called 'open' flows (see, e.g., Huerre and Monkewitz 
1990), for which the numerical instability simulations can be quite challenging because of 
the problem of boundary conditions on the inflow and outflow boundaries. 

It is worth remarking that the general forms of the basic flows, which have been obtained 
from the only requirement of separability of the corresponding stability problem, are reacher 
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than those remaining after specification to the exact solutions of the Navier-Stokes equa- 
tions. Thus, using the approach accepted in many stability studies, where the form of the 
basic flow is chosen quite freely to approximate the physical situation of interest, we could 
considerably enrich the list of relevant flows. However, our purpose was to provide exam- 
ples of a completely rigorous analysis that reduced the stability problem to an eigenvalue 
problem of ordinary differential equations. 

In addition, note that we have not yet exhausted the "direct approach" to separation of 
variables in the hydrodynamic stability equations. Changing the input data of the method 
given by equations ©-© may lead to new results. We considered here only the most natural 
generalization of the normal modes of the steady-state parallel flow analysis, which allows 
periodicity of perturbations in two new variables, and the order of the reduced equations 
was taken the same as that obtained in the parallel flow stability problem. If one uses 
other input data, it might be possible to construct new coordinate systems and basic flows 
providing separability of the system (|2j). 

We will also remark on the practical importance of specific basic flows that have been 
defined in the course of our analysis (some of them are discussed in the previous section). 
Those flows, mainly, are either ones over the stretching surfaces or the flows within porous 
channels possessing moving walls. The description of the flow near a stretching surface has 
many important applications in manufacturing processes in industry. A literature on the 
subject (see, e.g., the book by Pop and Ingham 2001) shows considerable research activities 
in this area. Solutions for physical situations, close to those considered in our Section III, 
can be found, for example, in Burde (1995a, 1995b), Youssef (1997), Nahapatra and Gupta 
(2003), Nazar et al. (2004). 

Laminar, incompressible and time-dependent flows that develop within a channel possess- 
ing permeable, moving walls have received considerable attention in the past due to their 
relevance in a number of engineering applications. Instances of direct application of such 
flows include the modeling of sweat cooling or heating, isotope separation, filtration, pa- 
per manufacturing, irrigation, and the grain regression during solid propellant combustion. 
From a different perspective, the sequences of expansions and contractions completed by 
channel walls enable a researcher to mimic more realistically peristaltic motion caused by 
pulsating walls and involving fluid absorbtion and filtration processes. For the cases, similar 
to those considered in the present paper, which pertain to a pipe that exhibits either injec- 
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tion or suction across porous boundaries while undergoing uniform expansion or contraction 
see, e.g., Uchida and Aoki (1977), Goto and Uchida (1990), Majdalani and Zhou (2003), 
Dauenhauer and Majdalani (2003). 

Thus, the exact solutions of the viscous incompressible Navier-Stokes equations deter- 
mined in this paper as the basic states, for which the linear stability problem is exactly 
separable, may be themselves of interest from both theoretical and engineering points of 
view. 
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